Search for evidence of two photon contribution in elastic electron proton data. 
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We reanalyze the most recent data on elastic electron proton scattering. We look for a deviation 
from linearity of the Rosenbluth fit to the differential cross section, which would be the signature 
of the presence of two photon exchange. The two photon contribution is parametrized by a one 
parameter formula, based on symmetry arguments. The present data do not show evidence for such 
deviation. 



I/-) i I. INTRODUCTION 

; 

Form factors (FFs) characterize the internal structure of composite particles. Since the first measurements 0, 
electromagnetic probes are traditionally preferred to the hadronic beams, because the electromagnetic leptonic inter- 
action is exactly calculable in QED, and one can safely extract the dynamical information about the hadronic vertex. 

1 ^ ■ Unfortunately, it is not a direct procedure, because, one has to introduce radiative corrections, which may become 

very large for the unpolarized cross section at large momentum transfer. Radiative corrections for elastic electron- 
hadron scattering were firstly calculated by Schwinger Q and are important for any discussion of the experimental 
determination of the differential cross section. They are also calculable in QED, within some assumptions concerning 

^vq , the hadronic interaction. There are standardprocedures applied to elastic eiV scattering data, but question may arise 
if the necessary accuracy has been reached j3| ■ 

f — , The largest value of momentum transfer squared, Q 2 , where the electric and magnetic proton FFs Ge p and Gmp 
have been extracted from the unpolarized cross section by the Rosenbluth method Q is Q 2 ~ 8.9 GeV 2 [j| and further 
extraction up to Q 2 ~ 31 GeV 2 assumes Ge p = Gmp/^p ((J>p = 2.79 is the magnetic moment of the proton). In 
order to show the sensitivity of the measurements with respect to Q 2 , FFs are often normalized to a dipole function: 



o 



G D (Q 2 ) - (1 + Q 2 [GeV 2 ]/0.7ir 2 . (1) 



i-S^ ' Recently, new developments, due to the very precise and surprising data obtained at the Jefferson Laboratory (JLab), 
in e+p — > e+p elastic scattering |7J,|8|, based on the polarization transfer method show that the electric and magnetic 
O ' distributions in the proton are different. 

The application of the polarization transfer method, proposed about 30 years ago 0, has been possible only recently, 
as it needs high intensity polarized beams, large solid angle spectrometers and advanced techniques of polarimetry in 
the GeV range. Experiments have been performed at JLab up to Q 2 = 5.6 GeV 2 and an extension up to 9 GeV 2 is 
• ^ in preparation 

The following parametrization for the ratio R of the electric and magnetic form factors well describes these experi- 
mental data 01: 



R = t*pG E p/G Mp = 1 - 0.13(Q 2 [GeV 2 ] - 0.04) (2) 

which implies that the ratio monotonically decreases and deviates from unity, as Q 2 increasing, reaching a value of ~ 
0.3 at Q 2 ~ 5.5 GeV 2 . 

Therefore, a clear discrepancy appears between the <5 2 -dependence of the ratio R of the electric to the magnetic 
proton form factors, whether derived with the standard Rosenbluth separation or with the polarization transfer 
method. This statement is confirmed by a reanalysis of the existing data [T^ and by recent measurements fl3l Il4|. 
This discrepancy is very puzzling, as no evident experimental bias has been found in the data or in the method used, 
and it has been sources of different speculations. 

One has to stress, at this point, that the observables are, on one side the differential cross section and on the other 
side, the polarization transferred to the scattered proton. The discrepancy is not at the level of the observables: FFs 
extracted with the polarization method, are not incompatible with the measured cross section: it has been shown that, 
calcualting Gmp from the measured cross section with the constraint on the ratio R from polarization measurements, 
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leads to a renormalization of Gmp of the order of 2-3% only, with respect to the Rosenbluth data [Tlj , well inside the 
error bars. 

Instead, the data are not consistent, if one refers to the slope of the e dependence of the reduced cross section, (e 
is the polarization of the virtual photon, e = [1 + 2(1 + r) tan 2 (0/2)] _1 ,O < e < 1 ) which is directly related to Ge p - 
The difference of such slope especially appears with respect to the last, precise data |14|. 

A possible question arises on the validity of the one-photon mechanism at large Q 2 , and, generally, on the radiative 
corrections to the differential cross section and to polarization observables in elastic eiV-scattering. If these corrections 
are large (in absolute value) for the differential cross section |15| . a simplified estimation of the radiative corrections to 
polarization phenomena [l6j shows that they are small for the ratio Pl /Pt of longitudinal to transverse polarization 
of the proton emitted in the elastic collision of longitudinally polarized electrons with an unpolarized proton target. 

In the standard calculations of the radiative corrections [l5j , the two-photon exchange mechanism is only partially 
taken into account considering the special part of the complicated loop integral, where one virtual photon carries all the 
momentum transfer and the second virtual photon is almost real. This contribution allows to overcome the problem 
of the 'infrared' divergence. But it has been pointed out 0] that, at large momentum transfer, the role of another 
part of the integral, where the momentum transfer is shared between the two photons, can be relatively increased, 
due to the steep decreasing of the electromagnetic form factors with Q 2 . This effect can eventually become so large 
(especially at large Q 2 ) that the traditional description of the electron- hadron interaction in terms of electromagnetic 
currents (and electromagnetic form factors) can become incorrect. 

Numerous tests of the possible 27 contribution for elastic ep scattering have been done in the past, using different 
methods: test of the linearity of the Rosenbluth formula for the differential cross section, comparison of the e + p and 
e~p-cross sections, attempts to measure various T-odd polarization observables, but no effect was visible beyond the 
precision of the experimental data. Only recently the non-zero T-odd asymmetry in the scattering of transversally 
polarized electrons, + p — > e~ + p has been detected [l8Hl9j . 

Following the arguments of Ref. |l7|. the relative contribution of two-photon exchange, compared to the main (i.e., 
one-photon) contribution should be visible at smaller Q 2 for heavier targets: d, 3 He, 4 He, because the corresponding 
form factors decrease faster in comparison with protons as Q 2 increases. One would expect to observe the two-photon 
contribution in ed-scattering at Q 2 > 1 GeV 2 , and for e TV-scattering at larger momentum transfer, Q 2 ~ 10 GeV 2 . 
In Ref. [2(j the possible effects of 27-exchange have been estimated from the precise data on the structure function 
A(Q 2 ), obtained at JLab in electron deuteron elastic scattering, up to Q 2 — 6 GeV 2 [2H |22|. The possibility of a 
27-contribution has not been excluded by this analysis, starting from Q 2 = 1 GeV 2 , and the necessity of dedicated 
experiments was pointed out. 

The complete calculation of the 27-contribution to the amplitude of the e^p — > e ± p-process requires the knowledge 
of the matrix element for the double virtual Compton scattering, 7* + N — > 7* + N, in a large kincmatical region 
of colliding energy and virtuality of both photons, and can not be done in a model independent form. However, one 
can partially reconcile the data, assuming that the radiative corrections contribute to the cross section linearly in 
e [23, 13 , or in framework of other model-dependent assumptions [2j| . These attemps give quantitatively different 
results, but have an effect in decreasing the slope of the differential cross section as a function of e. 

We previuosly showed [20I I2H I2H |28| . in a model independent way, that the presence of two-photon exchange 
destroys the linearity of the Rosenbluth fit, inducing a specific dependence of the differential cross section on the 
variable e, which is especially large for e — > 1. 

The purpose of this paper is to re-analyse the data about the differential ep scattering, with respect to the deviation 
of linearity of the Rosenbluth fit. The two-photon contribution is parametrized according to a simple one parameter 
form, which satisfies the necessary symmetry properties of the e dependence. A factorisation is assumed for the Q 2 
and e dependent terms. 



II. PARAMETRIZATION OF THE 2 7 CONTRIBUTION 

The traditional way to measure electromagnetic proton form factors consists in the measurement of the e dependence 
of the reduced elastic differential cross section a re d, at fixed Q 2 . Assuming that the interaction occurs through the 
exchange of one photon, one can write: 

« red (Q 2 , e) = e(l + r)[l + (225/m) sm 2 (8/2)} — ^- = rG 2 Mp (Q 2 ) + eG| p (Q 2 ), (3) 

OMott ail 

with r = Q 2 /(4m 2 ), m is the proton mass, E and are the incident electron energy and the scattering angle of the 
outgoing electron, respectively. 

The Rosenbluth separation Q allows to extract the electric and the magnetic form factors, from a linear fit to the 
reduced cross section, as a function of e, at a fixed Q 2 . The slope is directly related to Ge p and the intercept to Gmp- 
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However, due to the coefficient r, the weight of the magnetic contribution becomes larger, as the momentum transfer 
increases, reducing the sensitivity of the elastic cross section to the electric contribution. As an example, at Q 2 ~ 4 
GeV 2 , the term related to the electric FF contributes for less than 10% to the reduced cross section, assuming the 
dipolc scaling, whereas it would be as low as 2%, if one assumes the Q 2 dependence of R from Eq. One can 

estimate the level of precision required to extract the electric FFs from the experimental cross section at large Q 2 . 
In presence of 2j exchange, Eq. J3J) can be rewritten in the following general form: 

<Jred(Q 2 , e) = eG 2 E (Q 2 ) + tG 2 m (Q 2 ) + aF(Q 2 , e), (4) 

where a — e 2 / (47r), and F(Q 2 , e) is a real function (of both independent variables Q 2 and e), describing the effects of 
the ±7 (g) 27 interference. 

Note that, in the general case, F(Q 2 , e) contains two different contributions, one due to the effect of the two-photon 
contribution on the form factors Ge.m, and another one due to the third spin structure in the matrix element of 
eN scattering, induced by 27 exchange. Both contributions can be calculated only in framework of some model, 
considering different intermediate states for the 27 box diagrams. 

In such situation, any model independent statement concerning the function F(Q 2 ,e) is important, in particular 
concerning its e dependence, due to the large sensitivity of the extraction of G 2 E to the additional contribution F(Q 2 , e), 
at large Q 2 . 

We proved earlier [2(|, that the C-invariance and the crossing symmetry of hadron electromagnetic interaction, 
result in the following symmetry properties of F(Q 2 , e) with respect to the variable x - 



1 -e 

F(Q*,x) = -F{CF,-x). (5) 

This means that the I7 (^) 27 contribution F(Q 2 , e) has an essential non linear e dependence, at any Q 2 . For example, 
the additional (third) spin structure is generating the following e-dependence: 



F{Q\x) -+ e ^±£/^(Q 2 ,e) orF(Q 2 ,x) ^ i±lf^ (Q 2 , e ), (6) 

where the upper index, (T) or (A), corresponds to the tensor |27| or axial |2S| parametrization for the third amplitude. 

In order to estimate the upper limit for a possible 27 contribution to the differential cross section and the corre- 
sponding changing to G e ,m{Q 2 ), we analyzed four sets of data 0, 0, 0, 0], applying Eq. Q with the following 
parametrization for F(Q 2 , e): 



F(Q 2 ,e)^e ] /^f( T HQ 2 )- (7) 

It is important to stress that Eq. Q is a simple expression which contains the necessary symmetry properties of the 
I7 (^) 27 interference, through a specific (and non linear) e dependence. 
We checked that the parametrization: 



6 f( A )(n 2 \ 



F{Q\x)^ x j^-f^{Q 

gives qualitatively similar results. 

For the Q 2 dependence of /^ T, ' 4 - ) we take: 

r(T,A)tQ2\ _ ^ /g\ 

(1 + Q 2 [GeV] 2 /0-71) 2 (l + Q 2 [GeV] 2 /m 2 T A ) 2 ' 

where C is a fitting parameter, mx t A is the mass of a tensor or an axial meson with positive C-parity. Taking 
m 2 ^ A ~ 1.5 GeV 2 (typical value of the corresponding mass), one can predict that the relative role of the 27 contribution 
should increase with Q 2 . 

III. RESULTS AND DISCUSSION 

In presence of 27, the dependence of the reduced cross section on e can be parametrized as a function of three 
parameters, G E , G 2 M and C, according to Eqs. an( l ©• 
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Q 2 [GeV 2 



± A 



± A 



C±AC 



X 2 /n 



N ■ 

1 ' pot- 



Ref. 



1 


0.7 ± 0.4 


1.08 ± 0.04 


3 ± 3 




3 


f3] 


2 


1.1 ± 0.2 


1.03 ± 0.03 


1 ± 1 


0.1 


8 


f3] 


2.5 


0.7 ± 0.4 


1.06 ± 0.03 


2 ± 1 


0.1 


6 


m 


3 


1.1 ± 0.5 


1.02 ± 0.07 


1 ± 4 


0.1 


5 


m 


1.75 


1.0 ± 0.1 


1.05 ± 0.01 


-1 ± 2 


0.1 


4 


[5] 


2.5 


0.5 ± 0.3 


1.07 ± 0.02 


2 ± 3 


0.8 


7 


[5] 


3.25 


1.0 ± 0.3 


1.04 ± 0.03 


-3 ± 6 


0.1 


5 


u 


4.0 


0.8 ± 0.3 


1.04 ± 0.03 


1 ± 3 


0.1 


6 


[5] 


5.0 


0.9 ± 0.4 


1.02 ± 0.05 


± 5 


0.3 


5 


f5] 


0.65 


1.0 ± 0.2 


0.98 ± 0.03 


± 2 




3 


[131 


0.9 


1.1 ± 0.1 


0.99 ± 0.02 


-2 ± 2 




3 


[13] 


2.2 


1.2 ± 0.3 


1.03 ± 0.04 


-3 ± 3 




3 


[13] 


2.75 


0.8 ± 0.3 


1.06 ± 0.03 


0. ± 4 


- 


3 


[13J 


3.75 


1.6 ± 0.5 


1.0 ± 0.1 


-17 ± 15 




3 


[13] 


4.25 


1.8 ± 0.6 


1.0 ± 0.2 


-19 ± 26 




3 


[13] 


5.2 


4 ± 2 


0.8 ± 1.5 


-166 ± 195 




3 


[13J 


2.64 


0.9 ± 0.1 


1.05 ± 0.01 


± 2 


0.1 


5 


[14] 


3.2 


0.8 ± 0.2 


1.05 ± 0.02 


4 ± 4 


0.5 


4 


M 


4.1 


0.9 ± 0.4 


1.03 ± 0.05 


5 ± 12 




3 


[14j 



TABLE I: Form factors and 27-coefficient, Eq. 



In Fig. ^ we show a re d as a function of the variable e for different Q 2 values, together with the results of the three 
parameters fit, for the data from Ref. [5j. Fits of similar quality can be obtained for all the considered sets of data. 

In Fig. |2 from top to bottom, the electric and magnetic FFs, as well as the two photon parameter C, are shown 
as a function of Q 2 (solid symbols). The previously published data, derived from the traditional Rosenbluth fit are 
also shown (open symbols). The expectation from Eq. J2Jl is shown as a dashed line in Fig. Oltop). 

The numerical values of FFs and of the 27-coefficient for the data set analyzed here are reported in Table [I] The 
resulting parameter C is compatible with zero. The effect of including a third fitting parameter, the 27 term, is to 
increase the error on FFs. For the three points of Ref. at higher Q 2 , the coefficient C becomes quite large, but 
is still compatible with zero within the error. For this reason, these points are reported in the table but not in Fig. 
13 The present analysis does not give any evidence of a 27 contribution, as the possible 27 term always vanishes, 
within the error bars. The new values of Ge and Gm are compatible with the published values, deduced in frame of a 
standard e fit. The x 2 l n (where n is the number of degrees of freedom, i.e., the number of points TV minus the number 
of parameters which is equal to 3) is always smaller than unity. Varying A from 1 to 3 GeV 2 slightly modifies the 
coefficient C and its error, being always compatible with zero. In Fig. [21 (bottom) no systematic dependence of C 
with Q 2 is visible. The three points at the highest Q 2 from Ref. ^| are not shown in the figure due to their large 
error bar, which, in our opinion, prevents from any conclusion on a possible systematic trend. 

The average results on this parameter is < C >— 0.5±0.6. No deviation from zero is seen also for the individual 
sets of data, within the errors. Taking into account the a term which is included in Eq. (0J, the deviation of linearity, 
which is quantified by the parameter C is smaller than 1%. We also evaluated the relative weight of the different 
contributions to the cross section. We calculated the ratios of the magnetic, electric and two photon contributions 
to the reduced cross section and averaged for all the considered values of Q 2 and e. We find, as expected, that the 
dominating term is the magnetic term with a contribution at the level of (95.9 ±0.3)%. The electric term contributes 
for (3.8 ± 0.3)% and the average for the two photon term is (0.05 ± 0.05)%. 

Furthermore, no visible systematic effect appears, that could allow, at least, to guess the sign of the 27 coefficient. 
If one could determine at least the sign of the interference contribution, it would be possible to predict the relative 
value of the differential cross section for e~p and e + p-scattering. The following relation holds: 

a{e-p) - a(e+p) = <xF(Q 2 ,e) 

a(e-p) + a(e+p) eG 2 E (Q 2 ) + tG 2 m (Q 2 ) ' {) 
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Another approach - introducing ad — hoc a linear e contribution to the differential cross section - results in a corre- 
sponding decreasing of the G 2 E contribution, bringing the results in agreement with the polarization measurements. 
This can not be considered a proof of the presence of the 27 contribution, but shows that a small correction to the 
slope of the reduced cross section is sufficient to solve the discrepancy. Such correction may arise, for example, from 
a revision of the standard procedure of the calculation of radiative corrections, especially at e — > 1, which would 
produce a large effect on the value of G 2 E , as extracted from the Rosenbluth fit. 

The radiative corrections are taken into account in the data analysis, according to the prescription of Ref. 
Typically, the experimental results on the measured elastic cross section, are corrected by a global factor 5r: 

a = °Vnea S e - ^measi 1 + d R>- 

This factor contains a large e dependence, and a smooth Q 2 dependence, and it is common for the electric and 
magnetic part. At the largest Q 2 considered here, it can reach 30-40%. A linear e dependence can be considered a 
good approximation for the radiative corrections: 6r — eS' , but our conclusions hold for any monotonically increasing 
function. 

Therefore Eq. © can be rewritten as: 

<L{Q\ e) = ° re \Q\ - 6 R (Q 2 , e)} ~ G 2 M (Q 2 ) (r + ^^^y) [1 - ^'(Q 2 )} 



G 2 M {Q 2 



G 2 E (Q 2 ) 
G 2 M (Q 2 ) 



r8'{Q 2 ) 



(11) 



neglecting the term in e 2 . At a fixed Q 2 (dropping the Q 2 -dependence) we can see that the slope of the measured 
cross section has also a linear behavior with e, with a slope a: 



a = G%- tS'G 2 m . 



(12) 



This shows that the slope of the measured reduced cross section, as a function of e (which is related to the electric 



FF squared), can vanish, and, for 8' > 3 may even become negative. In Fig. [3] we show a rec i as a function of e, 

with and without radiative corrections, taken from Ref. One can see that the slope is compatible with zero, for 
the uncorrected data, starting from Q 2 < 2 GeV 2 , and, then, becomes negative. 

The extraction of Ge from the data requires a large precision in the calculation of radiative corrections, and in the 
procedure used to apply to the data. In particular at large e, an overestimate of the radiative corrections, can be 
source of a large change in the slope. 

A similar study has been done from the other sets of data and leads to similar results. 



IV. CONCLUSIONS 



From this analysis it appears that the available data on ep elastic scattering do not show any evidence of deviation 
from the linearity of the Rosenbluth fit, and hence of the presence of the two photon contribution, when parametrized 
according to Eq. (J7J). 

The new generation of experiments performed at large Q 2 makes use of large acceptance detectors and requires 
important corrections of the raw data for acceptance, efficiency, energy and angle calibrations 01 . Limits of the 
procedure and the approximations used for calculating radiative corrections were previously discussed in the literature 
I 

Radiative corrections at the large Q 2 reached in the considered experiments may reach 30-40%, whereas they are 
considered to be negligible in polarization eperiments (although one should stress that, in this case, no complete 
calculation is available). In particular they contain an important e dependence which has a large influence on the 
slope of the reduced cross section, changing even its sign, as shown in this work. The slope is directly related to the 
extracted electric form factor. 

A careful study of the effect of radiative corrections to the experimental data seems necessary, due to their large 
size and large e dependence, in particular for the large values of e, where the recent experiments have been performed. 
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FIG. 1: Reduced cross section at Q 2 =1.75, 2.5, 3.25, 4, and 5 GeV 2 . Data are from Ref. The lines are three-parameter 
fits according to Eqs. @ an d ©■ 
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FIG. 2: From top to bottom: Ge p /Gd, G m v / H P G d , and two-photon contribution, C according to Eq. Q with = 1.5 

GeV 2 . Data are from Refs. Q (squares); |f| (circles); jl3j (triangles); \\M (stars). The present results (including 27 contribution) 
are shown as solid symbols. The abscissa of the published data (open symbols) is shifted by -100 KeV, for clarity. The 
expectation from Eq. is shown as a dashed line. 
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FIG. 3: Reduced cross section with (solid circle) and without (open circles) radiative corrections, for Q 2 =1.75, 2.5, 3.25, 4, 5, 
6, and 7 GeV 2 . Data are from Ref. |5J. Two-parameters linear fits are shown as solid lines. The dashed lines show the slopes 
suggested by the polarization data. 



